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ABSTRACT 

We present the results of the VLT/VIMOS integral-field spectroscopic observations of the 
inner 28" x 28" (3.1 kpc x 3.1 kpc) of the interacting spiral NGC 5719, which is known 
to host two co-spatial counter-rotating stellar discs. At each position in the field of view, the 
observed galaxy spectrum is decomposed into the contributions of the spectra of two stellar 
and one ionised-gas components. We measure the kinematics and the line strengths of the 
Lick indices of the two stellar counter-rotating components. We model the data of each stellar 
component with single stellar population models that account for the afPe overabundance. We 
also derive the distribution and kinematics of the ionised-gas disc, that is associated with the 
younger, less rich in metals, more a-enhanced, and less luminous stellar component. They are 
both counter-rotating with respect the main stellar body of the galaxy. These findings prove 
the scenario where gas was accreted first by NGC 5719 onto a retrograde orbit from the large 
reservoir available in its neighbourhoods as the result of the interaction with its companion 
NGC 5713, and subsequently fuelled the in situ formation of the counter-rotating stellar disc. 

Key words: galaxies: individual (NGC 5719) - galaxies: kinematics and dynamics - galaxies: 
abundances - galaxies: spirals - galaxies: stellar content 



1 INTRODUCTION 

The presence of stars counter-rotating with respect to other stars 
and/or gas has been detected in several disc galaxies and is com- 
monly interpreted as the end result of a retrograde acquisitio n of 
extern al gas and subsequent star formation (see Bertola & Corsinil 
1 19991 for a review). Nevertheless, some special cases of counter- 
rotating stellar discs c ould have an internal o rigin induced by the 
presence of a bar (e.g.. lEvans & Colletllll994h . 

The demography of gaseous and stellar counter-rotating com- 
ponents in SO's and spirals is a key to understand their assem- 
bly process. The fraction of lenticular galaxies with a counter- 
rotating gaseous disc is consistent with the 50% that we ex- 
pect if all the gas in SO's is of external origin jSertoIa et al.l 
Il99^ . In contrast, less than 10% of them ho s t a d etectable 
fraction of counter-rotating stars jKuiiken et all Il996h . Large- 
scale counter-rotation is a rare phenomenon in spirals. In fact, 
less than 10% of the studied spiral galaxies host a counter- 
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rotating gaseous and /or stellar disc dKannappan & Fabricant|[200ll : 
IPizzella et al.l l2004h . The particular case of large-scale stellar 
count er-rotation in disc galaxies has been observed only in NGC 
4550 teubin et al. 1992; Rix et al. 1992; Emsellem etaL lOoj) 
NGC 7217 (Merrifield & Kuijken I9 94|), NGC 3593 Isertola et al.l 
1 1994 ICorsini et ah. 1998. ; .Garcia-Burmo et alj|2000l). N GC 4138 
( |joreetal.lll996h . and NGC 5719 jVergani et al.ll2007t). To inter- 



pret the observed frequencies of counter-rotations. IPizzella et al.l 
(2004) argue that the retrograde acquisition of small amounts of 
external gas can give rise to counter-rotating gaseous discs in gas- 
poor SO's only, while in gas-rich spirals the newly acquired gas 
is swept away by the pre-existing gas. Counter-rotating gaseous 
discs in spirals are formed only from the retrograde acquisition of 
amounts of gas larger than the pre-existing gas content. This is the 
case of the pur ely gaseous coun t er-rotating components det ected 
in NGC 3626 ('Ciri et al.' 'l995"; 'Garci'a-Burillo et al.' '1998") and 
NGC 4826 (Braun et al. 1994; Garci'a-Burillo et al. 2003). Counter- 
rotating stellar discs are produced by subsequent star formation. 
Therefore, in this scenario counter-rotating stellar discs are ex- 
pected to be made by younger stars than those of their host galaxy. 

This picture can be directly tested in th e NGC 5719/13 galax y 
pair, which has been recently studied by IVergani et all ( l2007h . 
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Figure 1. Fit of the galaxy spectrum (black) in the spatial bin at 10" East 
from the galaxy center. The best fit model {green) is the sum of the spec- 
tra of the ionised-gas component (magenta) and the two stellar components 
(blue and red). The latter are obtained convolving the synthetic templates 
vifith the best fitting Gaussian LOSVDs and multiplying them by the best 
fitting Legendre polynomial. The differences in the position of absorption 
line features and in the H/3 equivalent widths between the two stellar com- 
ponents (indicating different kinematics and stellar population content) are 
clearly evident. 



NGC 5719 is an almost edge-on Sab galax y with a prominen t 
skewed dust lane at a distance of 23.2 Mpc. Vergani et al ] l l2007h 
report a spectacular on-going interaction with its face-on She com- 
panion NGC 5713. Two H I tidal bridges loop around NGC 5719 
and connect it to NGC 5713 at a projected distance of 77 kpc. The 
neutral and ionised hydrogen in the disc of NGC 5719 are counter- 
rotating with respect to the main stellar disc. The kinematics of 
the ionized-gas disc and of both the counter-rotating stellar discs 
are measured out t o about 40" (4. 3 kpc) from the galaxy centre. 
In conclusion, Verg ani et al. propose a scenario where H I 

from the large reservoir available in the galactic surroundings was 
accreted by NGC 5719 onto a retrograde orbit and subsequently 
fuelled the in situ formation of the counter-rotating stellar disc. 

In this work, we will address the crucial piece of information 
which is still missing, i.e. proving that the stellar population of the 
counter-rotating disc of NGC 5719 is younger with respect to that 
of the stars in the galaxy main disc. 



2 OBSERVATIONS AND DATA REDUCTION 

The integral-field spectroscopic observations were carried out in 
service mode with the Very Large Telescope (VLT) at the European 
Southern Observatory (ESQ) in Paranal during dark time between 
28 April and 16 June 2009. The Unit Telescope 3 was equipped 
with the Visible Multi Object Spectrograph (VIMOS) in the Inte- 
gral Field Unit (IFU) configuration. The HR blue grism covering 
the spectral range 4150 - 6200 A and the 0'.'67 fibre-^ resolution 
were used. The instrumental spectral resolution measured at 5200 
A was 2.0 A (FWHM), equivalent to 115 km s^^. Observations 
were organised into three dithered on-target exposures of 2950 sec- 
onds each, alternated to three offset sky exposures of 280 seconds 
each. The average seeing measured by the ESO Differential Image 
Meteo Monitor was 1". 

Data reduction (bias subtraction, fibre identification and trac- 
ing, flat fielding, wavelength calibration and correction for instru- 
ment transmission) was performed using the VIMOS ESO pipeline 
version 2.2.1 (http://www.eso.org/sci/software/pipelines/). The dif- 
ferent relative transmission of the VIMOS quadrants was corrected 



by comparing the intensity of the night-sky emission lines. The 
three offset observations were used to construct three sky spectra. 
To compensate for the time variation of the relative intensity of the 
night-sky emission lines, we compared the fluxes of the sky lines 
measured in the offset and on-target exposures. The corrected sky 
spectra were then subtracted from the corresponding on-target ex- 
posures. Each exposure was organised in a data cube using the tabu- 
lated correspondence between each fibre and its position in the field 
of view. The three sky-subtracted data cubes were aligned using the 
bright galaxy nucleus as reference and co-added into a single data 
cube. 

In order to increase the signal to noise ratio (S/N), spec- 
tra from fibres mapping adjacent regions in the sky wer e added 
togeth er using the Voronoi binning method (Ca ppellari & CopinI 
l2003h . Some of the spatial bins were modified to include only 
spectra from the regions associated to intense H/? emission or the 
dust lanes crossing the galactic disc. This ensured us to minimise 
the contamination of the kinematic and stellar-population proper- 
ties to be measured from spectra obtained in regions where stellar 
counter-rotation was detected and chemical decoupling is expected 
jVergani et al]|2007h . In total, we have 130 spatial bins few arc- 
seconds square large. We tested the robustness of our results using 
different binning schemes. 



3 KINEMATICS AND LINE STRENGTH INDICES 

In order to measure the kinematics and stellar population proper- 
ties of the two counter-rotating stellar components in NGC 5719, 
we need to separate their contribution to the observed spectrum in 
each spatial bin, taking advantage of their different velocities that 
causes a wavelength shift of their spectr a. To do that , we modi- 
fied th e penalized pixel fitting code (pPX F,ICaDpellari"& Emselle^ 
l2004h in a similar way to that done by ' McDermid et al.l ( 1200^ . 
The code builds two synthetic templates (one for each stellar com- 
ponent) as linear combination of stellar spectra fro m the MILES 
librar y (at FWHM = 2.54 A spectral resolution, Beifiori et al.l 
I2OIOI) and convolves them with two Gaussian line-of-sight veloc- 
ity distributions (LOSVDs) with different radial velocities and the 
same velocity dispersion a. This as sumption is in agreement with 
the findings by I Vergani et al.l j2007h . Results are confirmed if two 
different velocity dispersions are considered, although the velocity 
fields and maps of stellar population properties appear more noisy 
due to the additional degree of freedom that increases the degener- 
acy between measured parameters. Gaussian functions are added to 
the convolved synthetic templates to account for ionised-gas emis- 
sion lines (H/3, [O in]AA4959, 5007, and [N I]A5198) and fit si- 
multaneously to the observed galaxy spectra. Multiplicative Leg- 
endre polynomials are included to match the shape of the galaxy 
continuum, and are set to be the same for the two synthetic tem- 
plates. Our techniqu e represents an improvement with respect to 
IVergani et al.l ( 120071) : in fitting the measured LOSVD with a dou- 
ble Gaussian, they assumed that the two counter-rotating compo- 
nents had the same stellar population. On the contrary, our syn- 
thetic templates account for different stellar populations since they 
are independently built from the MILES library. We check the 
robustness of our results using also the INDO-US Coude library 
of stellar spectra ( Valdes et al. 20 04) and the single stel lar popu- 
lation (SSP) synthesis models by IVazdekis et al.N2010h to build 
the synthetic templates. Figure [T] shows an example of the de- 
composition of the galaxy spectrum measured in one of the spa- 
tial bins where the two stellar counter-rotating components and 
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Figure 2. Velocity field of the main stellai' component (left panel), counter-rotating stellar component {central panel) and counter-rotating ionised gas (right 
panel) in NGC 5719. Centre is in R.A. = 14'M0™56.3=, Dec. = -00° 19'05.4" (J2000.0). Black contours: isophotes derived from the H/3 emission map. 
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Figure 3. Equivalent width of the Lick indices in the main (red diamonds) and counter-rotating (blue circles) stellar components. Predictions from single 
stellar population models bv lThomas et alj fcoiOi) are superimposed. Crosses indicate mean error bars associated to the equivalent widths. 



ionised gas are observed. We perform the double stellar compo- 
nent fit only in 73 spatial bins, where the separation of the two 
kinematic components is reliable, and do the single stellar compo- 
nent fit to the spectra of the other bins. The spectra of outermost 
spatial bins have a low 5'/A'', therefore only the ionised-gas emis- 
sion lines are fitted. The stellar kinematics are measured by tak- 
ing into account the difference in spectral resolution between the 
VIMOS spectra and the MILES library. We derive the properties 
of the stellar populations of the two stellar counter-rotating com- 
ponents of NGC 5719 by measuring the line stren gth of the Lick 
indice s H/?, Mgb, Fe5270, and Fe5335 as defined bv lWorthev et al.l 
( 1 19941) in the synthetic templates. To this aim, the spectra are set 
to rest-frame by adopting the measured radial velocity and con- 
volved with a Gaussian function to match the spectral resolu- 
tion of the Lick system (FWHM = 8.4 A). We also calculate 
the in ean iron index ( Fe) = (Fe5270 + Fe5335) /2 (Gorgas et al. 
Il99d) and the combined magnesium-iron index [MgFe] — 
\/Mg6 (0.82 • Fe5270 + 0.28 • Fe5335). The [MgFe]' index is al- 
most independent from g-enhancemen t and hence serves best as a 
metallicity tracer jThomas et ai]|2003h . 



We assume that the two synthetic templates are a good approx- 
imation of the mean stellar populations of the two counter-rotating 
components. This depends on the adopted library of stellar spec- 
tra and on the accuracy of the fitting program in recovering simul- 
taneously the kinematics (i.e, radial velocity and velocity disper- 
sion) and the population properties (i.e., the continuum and lines 
strength) of the two counter-rotating components. To test this, we 
carry out Monte Carlo simulations on a set of 1200 artificial galaxy 
spectra to test the reliability and accuracy of the procedure to mea- 
sure the kinematics and line strength of the Lick indices of the two 
counter-rotating components and quantify the errors. We generate 
each galaxy spectrum by summing two stellar spectra with differ- 
ent velocity separation (AV^), different fraction of stars in the main 
component with respect to the total (Fmain), different cr, and dif- 
ferent 5'/A'^. The two stellar spectra are chosen to have similar line 
strengths to those measured for the corresponding stellar compo- 
nent. Gas emission lines are also added to the artificial galaxy spec- 
tra. The modified pPXF routine is then applied to analyse the artifi- 
cial galaxy spectra as if they are real. Errors on the fitted parameters 
are estimated by comparing the input and output values. Errors are 
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Figure 4. Two-dimensional maps of age, [Z/H] and [a/Fe] of the main (upper panels) and counter-rotating (lower panels) stellar components, in the region 
where the two components are disentangled. Black contours: isophotes derived from the H/3 emission map. Centre is as in Fig.|2] 



assumed to be normally distributed, with mean and standard devia- 
tion corresponding to the systematic (As) and random (Ar) errors 
on the relevant parameters, respectively. We find that for typical 
values observed in NGC 5719 (AV > 150 km 30 < S/N < 
90, o- ~ 90 km s"\ 30% < i^main < 70%) the random errors 
on the recovered AV and Fmain are smaller than 25 km and 
11%, respectively, while systematic errors are negligible. The mea- 
surements of the line strengths are affected both by random and 
systematic errors: if the true equivalent width of the spectral index 
I of the two components are Ji and 1 2 with Ii < I2, the program 
measures (/i + As) ± Ar and (J2 — As) ± Ar. For the typical 
values observed in NGC 5719, random errors are 0.15 — 0.30 A, 
and systemic errors are 0.10 - 0.20 A, mildly depending on /i and 
l2- Figure [2] shows the measured two-dimensional velocity fields 
of the main and counter-rotating stellar components and ionised 
gas. Figure[3]shows the Lick indices [MgFe]', H/3, Mgb, and ( Fe) 
we measure in the bins where the two counter-rotating components 
are disent angled. The pr edictions of the single stellar population 
models bv lThomas eTal]f2010) are also shown for reference. Cor- 
rection for the offset to the Lick system is not applied, since no 
standard Lick stars were observed. Nevertheless, the VIMOS spec- 
tra are flux calibrated, therefore we expect the offset to be small 



and constant. This represents only a rigid shift to the points in Fig. 
[3] and therefore does not change our result. 

The luminosity-weighted metallicity ([Z/H]), Q-enhancement 
([a/Fe]) and population age are derived for each co mponent by a 
linear interpolation between the model predic tions by Thomas et al.l 
( I2OIOI) using an iterative procedure as done in llVIorelli et alj ( l2008h . 
We do not include in the interpolation all the spatial bins with 
equivalent width of the Lick indices located outside the model grid 
more than the mean errorbar. The two-dimensional maps of age, 
[Z/H] and [a/Fe] of the main and counter-rotating stellar compo- 
nents are shown in Figure |4] 



4 DISCUSSION AND CONCLUSIONS 

We present the results of the VLT/VIMOS integral-field spectro- 
scopic observations of the inner 28" x 28" (3.1kpc x 3.1kpc) of 
the spiral galaxy NGC 5719. The kinematics of the stars and ionised 
gas are very complex. The observed galaxy spectra are decom- 
posed into the contributions of three distinct kinematic components 
characterised by a regular disc-like rotation: one main and one sec- 
ondary stellar component and a ionised-gas component. The spec- 
tral decomposition is done using an implementation of the pPXF 
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routine, which allows to measure simultaneously the kinematics 
and stellar population properties of the two stellar components. 

The rotation of the main stellar component is receding to- 
wa rds East, like t hat of the stellar body of the galaxy as observed 
by IVerpani et alj i IzOOTI) out to ~ 40". We measure a maximum 
rotation velocity of ~ 150 km at about 10" (1.1 kpc) from 
the centre. The secondary stellar and ionised-gas components are 
counter-rotating with respect to the main stellar component, with a 
maximum rotation of ~ 200 km at about 10" from the cen- 
tre. We are able to resolve the two stellar components down to the 
innermost ^ 2" (0.2 kpc). At smaller radii, they have a too small 
velocity separation to be resolved. The median values of light frac- 
tion contributed by the two components in the spatial bins where 
they are resolved are Fmain = 56% and Fsecondary = 44% (with 
20% standard deviation), respectively. 

The ionised gas is detected all over the observed field of 
view. It is characterised by a strong H/3 emission, which is con- 
centrated in a twin-peaked morphology indicating an edge-on ring 
with a semi-major axis of ~ 7" (0.8 kpc), and an outer asymmet- 
ric m = 1 and fragmented spiral arc extending southwards with a 
semi-major axis of ~ 13" (1.5 kpc; see contours in Figs.|2]and|4l(. 
The two-stream fluid instab ility present in systems like NGC 3593 
jGarcia-Burillo et al.ll200(il) can favor this morphology. An hint of 
the inner ring is visible also in the H I posit ion-velocity diagram 
measured along the major axis of NGC 5719 jVergani et al.ll2007h . 
The secondary stellar component is associated to these features 
of the gas distribution, i.e. the counter-rotating stars are detected 
in a region enclosed by the ionised-gas structures. A similar phe- 
nomenon is observed in NGC 3593, where a concent rated gaseous 
ring is associated to the counter -rotating stellar disc jCorsini et al.l 
ll998l : lGarcfa-Burillo et al.ll2000h . 

The two counter-rotating stellar components are characterised 
by different chemical properties as it results from the measured line 
strength of the Lick indices (Fig. [3}. On average, the main stel- 
lar component has lower H/3 and higher [MgFe]', Mgi, and ( Fe) 
with respect to the secondary component. This translates imme- 
diately into different properties of their stellar populations (Fig. 
|4ll. The main stellar component has ages ranging from 2 to 13.5 
Gyr (median age = 4 Gyr with 4 Gyr standard deviation). It 
has nearly-solar metallicity (median [Z/H] = 0.08 dex with 0.20 
dex standard deviation). It displays super-solar enhancement (me- 
dian [a/Fe] = 0.10 dex with 0.10 dex standard deviation). The 
counter-rotating stellar population component is younger, with ages 
ranging from 0.7 Gyr to 2.0 Gyr (median age =1.3 Gyr with 0.6 
Gyr standard deviation). Its metallicity shows a radial gradient: it 
changes from sub-solar ([Z/H] ~ —1.0 dex) in the outskirts, to 
solar ([Z/H] = 0.0 dex) and super-solar ([Z/H] ~ 0.3 dex) in 
the centre. The youngest ages and highest metallicities are found 
in correspondence of the regions where the H/? emission is more 
intense. The Q-enhancement is super-solar (median [a/Fe] = 0.14 
dex with 0.05 dex standard deviation). 

These findings extend the results bv lNeff etal] j2005t) based 
on GALEX observations. They reported the presence of a young 
stellar component in the disc of NGC 5719 by analysing UV and 
optical images. Unfortunately, the poor GALEX angular resolution 
prevented them to associate it to the H/3 ring as we successfully do 
(Fig.|4}. Our estimate is in agreement with their lower limit to the 
actual age of such a young population (0.3 ± 0.1 Gyr). 

With our new observations and spectral decomposition tech- 
nique, we prove that the mean age of the counter-rotating disc, 
which is associated to the neutral and ionised gas disc, is indeed 
younger than the main stellar disc. This result shows that counter- 



rotating disc has been recently assembled. The median overabun- 
dance of the counter-rotating component (0.14 de x) indicates a 
star formation history with a time-scale of 2 Gyr (Thom as et al.l 
[2005) . More details about the assembly process could be derived 
by a further analysis of the east-west asymmetries measured in the 
maps of the stellar-p opulation properties. The scenario proposed by 
IVergani et al] j2007h that NGC 5719 hosts a counter-rotating stellar 
disc originated from the gas accreted during the ongoing merging 
with its companion NGC 5713, is finally confirmed. 
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